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Abstract 

In the large N c limit, the and A c can be treated as bound states of chiral solitons and 
mesons containing a heavy quark. We show that the soliton and heavy meson are bound 
in an attractive harmonic oscillator potential. The Isgur-Wise function for — ► A c e~ V e 
decay is computed in the large N c limit. Corrections to the form factor which depend on 
tjim/ttiq can be summed exactly (rriN and uiq are the nucleon and heavy quark masses). 
We find that this symmetry breaking correction at zero recoil is only 1%. 
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1. Introduction 



In the heavy quark limit, the form factors for semileptonic A b — > A c e v e decay 
are characterized by a single universal function r\(v ■ v') 



(A c (v\ a') | cY (1 - 75) b \A b (v, s)> = 77^ ■ «') ^'h" (1 ~ 7s) "(v, 



1.1) 



where t> M and -t/ M are the four- velocities of the and A c , respectively. The Isgur-Wise 
function rj(v ■ v') @ has logarithmic dependence on the heavy b and c quark masses which 
is calculable using perturbative QCD methods. The quark mass dependence can be put 
into a multiplicative factor 



rj(v ■ v') = C cb (v ■ v') rj (v ■ v'). 



;i.2) 



where 



and 



C cb (vv') 



~a s (m b ) 


-6/25 


~a s (m c )~ 


a L (vv') 


_a s (m c )_ 









a L (v ■ v') 
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[v ■ v r(v ■ v') — 1] 



r(v ■ v') 
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In ( v • v 



' + ^(v ■ v') 2 - l) . 



;i.3) 

;i-4) 

;i.5) 



y/(v ■ V') 2 - 1 

For very large heavy quark masses (and fi of order the QCD scale), C cb (v ■ v') has a 
rapid dependence on v ■ v' . The function i]q(v ■ v') is determined by low-momentum strong 
interaction physics. It depends on the subtraction point \x in a way that cancels the 
subtraction point dependence of C cb (v ■ v'). At zero recoil, i.e. v ■ v ' = 1, 770 is independent 
of fx and is normalized to unity @] @ @ by heavy quark flavor symmetry, 



Vo(l) = 1. 



;i.6) 



Some of the low momentum properties of QCD are determined by its symmetries (e.g. 
chiral symmetry and heavy quark symmetry). Those nonperturbative aspects of the theory 
which are not determined by symmetries cannot be treated using perturbation theory in 
the strong coupling constant. QCD, however, does have an expansion parameter which 
can be used to study low momentum features of the strong interactions analytically. In 
the limit that the number of colors N c is large, the theory simplifies and many predictions 
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are possible |J. The main purpose of this paper is to examine rjo(v ■ v') in the large N c 
limit. 

In the large iV c limit, baryons containing light u and d quarks can be viewed as soli- 
tons [Q of the nonlinear chiral Lagrangian for pion self-interactions. Baryons containing 
a single heavy c (or b) quark and light u and d quarks are then described as bound states of 
these solitons with D and D* mesons (or B and B mesons) [f)|[10]. The large iV c behavior 



of rjo(v ■ v') can be determined using the bound state wavefunctions of the and A c . 
For large iV c , the function r/o (f • v') is strongly peaked about zero recoil since any velocity 
change must be transferred to ~ iV c light quarks. Independent of the details of the bound 
state approach, we find that 

r] {v ■ v') = exp[-XN^ 2 (v ■ v' - 1)], (1.7) 

where A is a constant of order unity. This equation is valid for v ■ v — 1 of order iV c 
In this kinematic region, 770 falls from unity to a very small quantity. The derivation of 
this result is the main purpose of this paper. In addition, the effect on 770 of corrections 
to the heavy quark limit that depend on m^lmQ is examined. Some features of the 
soliton picture of heavy baryons not discussed in previous work on this subject will also be 
derived here. In Refs. [|IU| [ITJ it was shown that the leading term in the chiral Lagrangian 
for heavy-meson-pion interactions gives rise to a heavy-meson-soliton potential that is 
attractive at the origin in the Aq, Eq and Eg channels. We show in this paper that the 
curvature of the soliton-heavy-meson potential is positive, indicating that the origin is a 
stable minimum of the potential energy for these channels. The curvaturve is negative for 
the exotic channels. 



2. Aq as a Heavy Meson-Soliton Bound State 

The starting point for discussing soliton-heavy-meson bound states is the chiral La- 



grangian for the interactions of mesons containing a heavy quark Q with pions [|12]. In 
the limit itlq — > 00, the total angular momentum of the light degrees of freedom, Si, is a 
symmetry generator. The lowest mass mesons with Qq a (qi = u, qi = d) flavor quantum 
numbers have si = 1/2 and form a degenerate doublet consisting of pseudoscalar and vec- 
tor mesons. In the case Q = c, these are the D and D* mesons, and in the case Q = b, 
these are the B and B mesons. 
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It is convenient to combine the fields P a and P* M for the ground state se = 1/2 mesons 
into the bispinor matrix 



1 + 



[P*a^ - Pall] 



(2.1) 



where is the heavy quark four velocity, and v 2 = 1. The vector meson field is constrained 
to satisfy v^P* = 0. In this section, we work in the rest frame of the heavy meson, 



(1,0). Under the heavy quark spin symmetry, 



H a — > SH ai 



(2.2) 



where S G SU{2) V is the heavy quark spin transformation. The transformation property 
of H under SU(2)l x SU(2)r chiral symmetry has an arbitrariness associated with field 
redefinitions. We will use the basis chosen in Ref. [1T0|1* with the transformation rule 



(2.3) 



under SU(2)l x SU(2)r, where R G SU(2)r. It is also convenient to introduce the field 



The Goldstone bosons occur in the field 

(2iM\ 

where 

M -[ 7T- -7T°/V2_ 

and / ~ 132 MeV is the pion decay constant. Under SU(2)l x SU(2)r 



2.4) 



with L G Sf7(2) L and R G SU(2) R . Under parity, 



(2-5) 



(2.6) 



(2.7) 



E(x°,x) -> E f (x ,- 



(2- 



since M(x°,x) — > — M(x°, — x). If transforms under chiral symmetry as in Eq. (27 3), 
then the parity transform of iif must transform under chiral symmetry with a factor of L'. 



The computations are repeated for the £ basis in the appendix. 
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Consequently, in the basis we are using, the action of parity on the H field is somewhat 
unusual, 

H a (x°,x) ->7 fr 6 (x ,-5)7 E+ 6 a (x ,-5). (2.9) 
The chiral Lagrangian density for heavy meson-pion strong interactions is 

£ = -i TrHv ■ dH + - TrM/E f ^E + — TrHH^^d^ + . . . , (2.10) 

where the ellipsis denotes the contribution of terms containing more derivatives or factors 
of 1/rriQ. The coefficient of the second term in the Lagrangian density Eq. ( |2.10| ) is fixed 
(relative to the first) by parity invariance. The coupling g determines the D* — > Dtt decay 
rate. Present experimental information on the D* width and the D* — ► Dn branching 
ratio implies that g 2 < 0.4 |T3|]. The constituent quark model predicts that g is positive. 

The soliton solution of the 577(2) l x SU{2)r chiral Lagrangian for baryons containing 
u and d quarks is 

E = A(t) X (x)A-\t), (2.11) 

where 

E = exp (iF(r) x ■ f) , (2.12) 

and r = A(t) contains the dependence on the collective coordinates associated with 
rotations and isospin transformations of the soliton solution. For solitons with baryon 
number one, F(0) = —tt and F(oo) = 0. The detailed shape of F(r) depends on the chiral 
Lagrangian for pion self interactions including terms with more than two derivatives. We 
expect that T,q(x) has a power series expansion in x. Consequently, the even powers of 
r must vanish when F(r) is expanded in a power series in r, e.g. F"(0) = 0. The chiral 
Lagrangian for pion self interactions is of order N c . However, the chiral Lagrangian for 
heavy-meson-pion interactions is only of order one. Thus, to leading order in iV c the shape 
of the soliton F(r) is not altered by the presence of the heavy meson. 

In the large N c limit, baryons containing light u and d quarks are very heavy and time 
derivatives on the S field can be neglected. Consequently, it is the interaction Hamiltonian 

H T = J d 3 x TrHH^^djZ + . . . , (2.13) 



with E given by Eqs. (p.ll[) and (|2.12j ) that determines the potential energy of a configu- 



ration with a heavy meson at the origin and a baryon at position x. Neglecting operators 
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with more than one derivative (the ellipsis in Eq. (j2.13|) ), and expanding the interaction 
potential operator in x gives 



Vi(x) = g S J £H I k H Tr At* A- 1 r k { 8 ij F'(0) - § r 2 [F'(0)f + ±r 2 F"'(0) 



+x l x J 



[F'(0)f + \F'"(0) + e ijm x m [F'(0)} 2 \ + 0(x 



(2.14) 



where Sj H denotes the angular momentum of the light degrees of freedom of the heavy 
meson, and 1^ denotes the isospin of the heavy meson. The interaction potential has terms 
which superficially have the wrong parity, e.g. the term involving the e symbol. However, 
these terms are required because of the £t factor in the parity transformation of H in 
Eq. 0. 

The Aq baryon has isospin zero and total angular momentum of the light degrees 
of freedom equal to zero. In the large iV c limit, it arises from a bound state of nucleons 
with P and P* mesons. Baryons with I > 1/2 such as the A cannot produce a heavy 
baryon bound state with 1 = 0. On nucleon states, Tr AT % A~ X T k is equal to — 8S l N I%/2> 
where S l N is the spin of the nucleon, and I N is the isospin of the nucleon [[nj. Using this 
simplification, the potential operator becomes 



Vt(x) = Vj {0} + V I (L > + V^> + 0(x 



(D T>(2) 



(2.15) 



(n) 

where Vj denotes the term of order r n in the potential 



V, 



(o) 



■UF'(0)I H -I N S eH -S 



V} (1) = § g [F'(0)Y Ih-Inx- (S m X S N ), 



v T (2) = -% g r H .r N {Sin ■ s N \-y [F\o)r + yF">(o) 

+ (S lH ■ x)(S N ■ x) [f [F'(0)] 3 + |F W (0)] } . 



(2.16) 



The potential V/(af) commutes with the total angular momentum of the light degrees of 
freedom, L + Sm + Sn, where L is the orbital angular momentum. 

To find the Aq wavefunction ^\(x) and its potential energy V\(x), the potential 
energy operator must be diagonalized at each point x on the product space of nucleon 
heavy meson states (e.g. \p, f) \P* , ~ ))• It is convenient to consider linear combinations 
of these product states that have definite isospin I = Ih + In, spin of the light degrees of 
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freedom, Si = Sm + Sn, an d total spin S = Sq + Si. These states are labeled \I, s, si). 
It is straightforward to diagonalize Vi(x) in this basis, and we find that 



1-F'(0)x-(Si H xS n ) |0,i,0> 



x 



[2.17) 



is an eigenstate of Vi(x) with eigenvalue 



V A (x) 



-lgF'(0)+gr 2 [i [F' 
■%gF'(0) + Ur 2 , 



^'"(0) 



where k is defined by 



« = 9 



| [F'(0)f 



5 TTlW 

6^ 



(0) 



(2.18) 



(2.19) 



The form of the wavefunction in Eq. (|2.17|) can also be found using Eq. (|2.9| ) and demanding 
that it has definite parity. The factor in square brackets in Eq. fl2.17f) compensates for the 
factor of £t j n Eq. flp). In Eqs. (|2T7|) and flTl8|) , rA QCD is treated as a small quantity* 
and the wavefunction is given to linear order in tAqcd? while the potential energy is 
given to quadratic order in tAqcd- The spatial part of the wavefunction <p(x) has a more 
rapid dependence on r which will be computed later in this article. As r goes from zero to 
infinity, F(r) goes from — tv to 0. Consequently, we expect that F'(0) is positive and F"'(0) 
is negative. This is true for example for the solution given in Ref. [|14| where a particular 
four derivative term in the chiral Lagrangian for pion self-interactions is used to stabilize 
the soliton. Furthermore, the constituent quark model suggests that g is positive. Thus, 
Eq. ( ]2.18p implies that the Aq is bound by a harmonic oscillator potential with k > 0. 

In the limit N c —>■ oo, the nucleon is infinitely heavy and terms in the Hamiltonian in- 
volving the nucleon momentum are neglected. The lowest energy state then has the spatial 
wavefunction <p(x) = 5 3 (x) corresponding to the minimum energy classical configuration 
where the nucleon is located at x = 0. The 1/N C terms in the Hamiltonian involving the 
nucleon momentum give the spatial wavefunction a finite spread and it is this finite extent 
which is responsible for the Isgur-Wise function rjo(v ■ v'). Terms of order 1/A^ C involv- 
ing the momentum are found by writing E = E(x — r(t)) and quantizing the collective 
coordinate r (t). This procedure yields 



#kin = 7T t F '(°) (Ih-In 

2m N 3 V 



S 



N -P 

m N 



(2.20) 



Aqcd denotes a nonperturbative strong interaction scale that is finite as _/V c — > oo. 
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where is the mass of the nucleon. The first term is the usual kinetic energy of the 
soliton and the second term in Eq. ( |2.20| ) results from the second term in Eq. ( |2.1U| ). The 
Schrodinger equation is 

[#kin + V A (x)\ # a (£) = E^f A (x), (2.21) 
which implies that (j)(x) obeys the differential equation 

V 2 



2m n 



+ V A (x) 



x)=E(f>(x), (2.22) 



in the rest frame of the bound state. Note that the term linear in p in Eq. ( |2.2Up is 
necessary for the wavefunction in Eq. (|2.17| ) to obey the Schrodinger equation (|2.21|) . This 
result is not surprising because both the term linear in p and the part of the wavefunction 
linear in x arise from the peculiar definition of parity in Eq. ( |2.9| ). In deriving Eqs. ( |2 . 2 1| ) 
and ( |2.22|) we have treated the term linear in p in H^ n as a perturbation and neglected 
its action on the "small" part of the wavefunction \I/a(x) (i.e. the piece proportional to 
x). As we shall see shortly, for large N c the term proportional to x in ^>a(x) an d the term 
linear in p in -ffkin are subdominant and can be neglected in the calculation of r/o (v ■ v'). 
In general, for large N c we expect the potential V\(x) to have the harmonic oscillator 

form 

V A (x) = V + \kx 2 . (2.23) 

The absence of a term linear in r requires F"(0) = which is a consequence of the derivative 
expansion of the chiral Lagrangian for pion self-interactions. Quantum corrections can 
induce nonanalytic behavior in x, but because of an explicit factor of 1/N C such terms are 
less important than those we have kept. The particular expression for k in Eq. (|2.19| ) is, 
however, model dependent and arises from keeping only terms with one derivative in the 
chiral Lagrangian for heavy meson-pion interactions. 

The model independence of the results of this paper follows from an analysis of large 
N c power counting. In the limit uiq — > oo, the typical size (or momentum) of the bound 
state wavefunction occurs when the kinetic energy p 2 /2m n and potential energy kx 2 /2 
of the bound state contribute equally to the total energy E — Vq, 

t ~ (ktyi tv) , p ~ (nrriN) 1 ^ 4 ■ (2.24) 

Since k is of order Aq CD and is of order Aqcd-^ c , the A^ c dependence of the typical 
binding energy is given by 

E — Vq ~ Aqcd-^T 1 / 2 . (2.25) 
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It is now straightforward to see that higher order terms in the effective Lagrangian can 
be neglected. Any term in the effective Lagrangian is a function of r = tn — th, nv, 
and fjj- Terms involving only the soliton field carry an overall factor of iV c . These terms 
are independent of the relative coordinate r and depend only on rjv- In the effective 
Hamiltonian, the dependence of these terms on enters through powers of the nucleon 
momentum. A term with n powers of the nucleon momentum has the following large N c 
behavior, 

AqcdNc {w^) ~ A QCDN c N~ 3n /\ (2.26) 

Thus, all terms with higher powers of the nucleon momentum than the leading order kinetic 
energy term p 2 jlmjq are suppressed by more powers of 1/N C and can be neglected in the 
large N c limit. Terms in the Lagrangian which involve the interaction between the soliton 
and the heavy meson can depend on r, but they are at most of order one in the large N c 
limit. The typical scale of the r dependent interaction is Aqcd 5 so higher order interaction 
terms have the form 

Aqcd (A QCD r) m (j^J " „ AgcD^X -8 " 74 - (2.27) 

Hence, all the higher order interaction terms in the effective Lagrangian other than the 
harmonic potential are higher order in 1/N C and can be neglected (including the term 
linear in p in if kin Eq. (|2.2C| )). 

There is also a class of 1/mg corrections which can be summed exactly. So far, we 
have concentrated on the order of limits uiq — ► oo followed by iV c ^ oo. We now switch 
to the situation in which rriQ,N c — > oo simultaneously, with the ratio Aqcd A^/mQ held 
fixed. The only additional term in the effective Lagrangian which must be included in this 
double scaling limit is the kinetic energy of the heavy meson, 

p 2 _ m N p 2 _ ( AqcdA^ p 2 ^ 22g ^ 



m H niQ m N \ uiq J m N ' 



which is of the same order as a term we have included, since Aqcd N c /rriQ is of order one. 
Higher order terms in 1/mg such as 



4 / \ 3 4 

V I m N \ V 



m Q \ m Q J m N ' 



(2.29) 



are of order one times terms which can be neglected by the power counting arguments of 
the previous paragraph, so they too can be neglected. Additional 1/mg effects arise from 
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1/toq operators in the heavy quark effective theory. The P* — P mass difference is of 
order Aq CD /trq ~ (A^ c Aqcd/wq) Aqcd/A^ which is subleading in N c . Higher derivative 
operators in the current are also suppressed. For example, the leading correction 

— Tr H (c) ^ (1 " 75) ~ AQCDiVcl/4 ~ f at-3/4 5 (2 .3o) 

since the typical momentum of the heavy quark in the baryon is of order Aqcd-^c^ 4 - Thus, 
no additional terms other than the kinetic energy of the heavy meson are relevant. 



3. Ab — ► A c e v e Decay 

For non-relativistic A c recoil, the matrix element of the weak current Eq. ( |1 . 1| ) in the 
A?, rest frame is 

(A c (v',s')\c^(l- l5 )b\A b (v,s)) = j d 3 p' I d 3 p <P:(P')MP) 
1 (iv(-^ + m N v' , s') I N(-0, s)} (D(f + m D v' )\ (1 - 75) b \~B[p)) + ... 

where the ellipsis represents terms involving at least one vector meson. <p c and (pb are the 
Fourier transform of the ground state wave function 

^{p) = ^2^^)3/8 ex P (-P 2 /2VM^«) > ( 3 - 2 ) 



where k is defined in Eq. (|2.19|) , and \i c ^ is the reduced mass \i = run / [t^n + ^h] 



of the bound state, with m# = m^, for the c and b subscripts, respectively. The Aq 
state is a superposition of products of N with P and P* states. However, we do not need 
the details of the Clebsch-Gordan structure* of the Aq state, as will become clear soon. 
The nucleon matrix element vanishes unless 

V = p + m N v\ (3.3) 



at which point, for the term explicitly displayed in Eq. (|3.1|) , the required heavy meson 
matrix element is 

(D (p+ (m N + m D )v' )|c 7 M (l -75)6 \B (p)> . (3.4) 



* The factor of 1/4 is the square of the appropriate Clebsch-Gordan coefficient for the term 
displayed explicitly in Eq. fl3.ip . 
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This matrix element can be evaluated in terms of heavy meson form factors. In the Aq 
heavy meson-nucleon bound state, the typical momentum is of order (p) ~ {m^K) 1 ^ 
and so form factors for semileptonic Aj, — > A c e~ V e decay are smooth functions of 
ttinv'/ (tunh) 1 ^ 4 which is of order Nc v' ■ Consequently, we are interested in the kine- 
matic region v of order N c , and so for large N c we can replace the heavy meson form 
factors by their rapidly varying part, C c b(v ■ v'), and neglect the slow variation in the 
meson Isgur-Wise function. This is why we do not need the details of the Clebsch-Gordan 
structure of the Aq bound state. The A5 — > A c e~ V e form factor in the large N c limit can 
be written in terms of an Isgur-Wise function, as in Eq. ( |1 . 1|) , with 



no 



d 3 p (p* c (p + m N v') 4> b (p 



1/4 



2 (HcLib) 

V^b + V^c 



3/2 



exp [-m 2 N v' 2 /2y/K (yfji b + y/Ji c )] 



(3.5) 



Since for non-relativistic recoils v' 2 ~ 2(i> -v' — 1), the expression for r]o in a general frame 
has the form 



Vo(v ■ v') 



2 (pel**) 
V^b + V^c 



i/4-l 3 /2 



exp [-m 2 N {v-v' - l)/y/K(y/jl b + y/fi c )] ■ (3.6) 



The result we have obtained is valid in the heavy quark and large N c limits, where we 
have included all terms of order (mjv/mQ)™ ~ (A? c Aqcd / 'wg) n • The baryon form factors 
for Ab — > A c e~V e are still written in terms of a single function in this limit, even though 
we have included a class of 1/mg corrections. In the large A^ c limit, the function tjq(v ■ v') 
falls off rapidly away from zero recoil. Derivatives of 770 at zero recoil diverge as N c — ► 00. 
Eq. ( |3.6|) indicates that the m th derivative is of order N 3171 ^ 2 , and includes all contributions 
of this order neglecting less divergent pieces. For example, there could be corrections to 

1/2/ 

rjo of the form N c (y -v — 1) in the exponent. This term is small compared with the 



3 /2 

leading term N c (y • v — 1), but is significant when v ■ v — 1 is of order A^ c 



-1/2 



Eq. ( |3.6| ) is valid in the large N c limit in the region where v ■ v' — 1 < 0(N C 3 ^ 2 ) 
the region where r/o falls from about unity to a very small quantity. 
At zero recoil in the limit that to& —>■ 00, 770 has the form 



Thus 
i.e. in 



7?o(l) 



2 {^ c m N 



,V4 



1 3/2 



N 



1 - -2- 
X 64 



m N 
m D 



+ ... 



(3.7) 
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where we have expanded the exact expression in a power series in m^/m£,. The term 
linear in ra^jrar, vanishes, which is consistent with Luke's theorem [|15]. For physical 
values of /rriD, the correction to the symmetry limit prediction 770 (1) = 1 is only 1%. 

The parameter k can be determined to be (530 MeV) 3 in the Skyrme model using 
the shape function used in Ref. ||14j| , and the value of g obtained in Ref. ||11|| . With this 
value of we find that the orbitally excited Aq state should be about 400 MeV above the 
ground state, and that the form factor 770 of Eq. ( |3.6|) is 



r] (v ■ v') ~ 0.99 exp [-1.3 (v ■ v' - 1)] , (3.8) 

using the known values of m^v, mp and mg. The Skyrme model prediction for k is 
sensitive to the precise shape of the soliton solution. A better way to determine k is to 
use the experimentally measured excitation energy of the orbitally excited Aq, which is 

The large N c predictions of this paper rely on the number of light quarks in the heavy 
baryon being large. For N c = 3 there are only two light quarks in Aq, so we expect our 
results to only be qualitatively correct. Nevertheless, it is interesting that the baryon form 
factors are calculable in the large N c limit of QCD. There are other results that can be 
computed using the methods developed here. For example, the Isgur-Wise functions for 
transitions to excited states are also calculable. It should also be possible to derive our 



results using the methods of Witten [16]. It would be interesting to explore that approach. 
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Appendix A. The £ Basis 

In this appendix, we briefly discuss the computation of the binding potential in the £ 
basis discussed in Ref. [[TOf. The notation is the same as found in Ref. [|10|. The £ basis is 
singular at the origin, but has a simple transformation rule for the H field under parity, 

H(x°,x) -> <y°H(x°, -x)7°. 
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The interaction Hamiltonian in this basis is 

Hj = -| J d 3 iTVMf 7 5 (tfdjt - £3^) . (A.l) 

Expanding the Goldstone boson field 

£ = exp(iF(r)x-f/2), (A.2) 

in a power series about x = 0, and using £ = A£ we get the interaction potential 

V7(x) = 2s Tr Ar^r* |f'(0) ■ x I k H & - \Sl H I k H ) 

+r 2 [*"(0)f - &*""( )) S} H 4 + &ff " x 4 (-i [F'(0)] 3 + §F"'(0)) } . 

(A.3) 

Note that only even powers of x occur in Eq. (|A.3| ) because the parity transformation 
is trivial in the £ basis. To leading order in N c , the kinetic term of the soliton can be 
neglected, so that states with a definite value of A are eigenstates of the Hamiltonian. It 
is convenient to consider the soliton state that is an eigenstate of A with A = 1. States 
with other values of A can be obtained by an isospin transformation, and so have the same 
energy. On states with A = 1, the interaction potential reduces to 

V T (x) = 4g (F'(0) (S tH -x Ih-x- \S m ■ I H ) 

+r 2 (i [F'(0)f - iF"'(0)) Sih -I H + S m -xI H -x(-± [F'(0)f + §F'"(0)) } 

(A.4) 

This Hamiltonian is singular at the origin because of the coordinate singularity in the £ 
basis, and the eigenstates of the Hamiltonian will also be singular at the origin. We know 
that the singularity at the origin has the form f ■ x, because that is the transformation 
function from the singular £ basis to the non-singular basis used earlier in this article. We 
therefore write the eigenstates of the interaction potential Eq. ( |A.4| ), in terms of new 
eigenstates \<j>) which are related by the unitary transformation 

M =<2I H - x) |0) . (A.5) 

The interaction potential in the \4>) basis is 

Vi(x) = (2 I H ■ x) V T (x) (2 I H ■ x) 

= 2g {F'(0) S m ■ I H + r 2 (-§ [F'(0)f + ^F"'(0)) S m ■ T H (A.6) 
+ (J [^'(0)f + §^"'(0)) [S IH -xIh-x- \r 2 S lH ■ I H ) } , 
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where we have used the relations 



Th-x) 2 = \ 



4' 



(A.7) 



(f H ■ xj (Si> H ■ I H ^j (T H ■ = \ [Sm ■ (Jh ■ xj - \ (§iH ■ Ih 
that follow from the anticommutation relation 

for the isospin-1/2 operators Ih- It is convenient to define the operator Kh = Ih + SgH, 
in terms of which Eq. ( A.6 ) can be rewritten as 

Vi(x) = g [F'(0) (k h - §) + r 2 (-± [F'(0)] 3 + 4^(0)) (it* - f 
+ (i [F'(0)] 3 + ±F"'(0)) (Kh-xKh-x- \r 2 K H ) } . 



(A.9) 



The allowed values of Kh for the H field are Kh = and Kh = 1 obtained by combining 
Ih = 1/2 with Sm = 1/2. The states with Kh = are the bound physical baryon states 
containing a heavy quark. On the Kh = states, Eq. (|A.9|) reduces to the potential 



Vf = -%g{F'(0) + r 2 (-1 [F'(0)] 3 + AF"'(0))} 

= -|^F / (0) + i K r 2 , 

where k is defined in Eq. ( |2.19| ). Thus the states with Kh = are bound in an attractive 
harmonic oscillator potential, and the potential at the origin has the value — 3gF' (0)/2. 
The potential is more complicated for exotic states which have Kh = 1- The last term in 
Eq. ( |A.9j ) is an irreducible tensor operator with Kh = 2 and contributes a spin-orbit term 
to the interaction potential (but note that it does not cause any mixing between Kh = 
and Kh = 1 states). The spherically averaged potential for Kh = 1 states is simple to 
compute, 

Vj Kh=1 (x) = y{F'(0) +r 2 (-i [F'(0)] 3 + ±F»'(0i 



lg F >(0) - l« r 2 



(A.11) 



2^ 6 

Thus the Kh = 1 states are unbound since the potential at the origin is positive, and the 
interaction potential is a repulsive inverted harmonic oscillator potential. The Aq state is 
obtained from the Kh = state by applying a projection operator [11]] , so the interaction 
potential for the Aq state is Eq. ( |A.10| ) . This is precisely the potential given in Eq. ( |2.18| ) 
of the text. 
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The kinetic term in the £ basis on the A = 1 soliton states has the form 

1 -2 2 • / -» \ 
£ = -Mjvf - —x- [x x 7 H J . (A.12) 

The first term is the usual soliton kinetic energy, and the second term is from the expansion 
of 

C = ~ TrW (^d^ + ea M e f ) • (A.13) 
The kinetic Hamiltonian in the £ basis obtained from the Lagrangian Eq. ( |A.12[ ) is 

1 ( ^ 2 . -, \ 2 



^ kin= 2M^^ + ^ fX/ ^ • (A ' 14) 
Transforming from the singular basis using Eq. ( |A.5| ) gives the kinetic energy 

-»2 

Hi in ={2I H . x) H kin (2 T H -x)= JL-. (A.15) 

The bound state problem in the £ basis reduces to that of a three-dimensional harmonic 
oscillator with a conventional kinetic term. 
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